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1. Introduction 

Gravastars (gravitational-vacuum stars) are hypothetical objects mooted as 
alternatives to standard Schwarzschild black holes [H [2]. Typically the interior is 
some simple nonsingular spacetime geometry such as de Sitter space, the exterior 
is some close approximation to the Schwarzschild geometry, and there is some 
complicated transition layer near the location where the event horizon would 
otherwise have been expected to form. Thus, in the traditional gravastar picture, 
the transition layer replaces both the de Sitter and the Schwarzschild horizons, and 
consequently the gravastar model has no singularity at the origin and no event 
horizon, as its surface is located at a radius slightly greater than the Schwarzschild 
radius. In this model, the quantum vacuum undergoes a phase transition at or near 
the location where the event horizon is expected to form. Considerable attention 
has been devoted to these objects, and to closely related "dark stars" , "quasi black 
holes" , "monsters" , "black stars" , and the like O IH O [6] . Related models, analyzed 
in a different context, have also been considered by Dymnikova [7]. Some models use 
a continuous distribution of stress-energy, which on rather general grounds must be 
anisotropic in the transition layer [HI [9] . Other models idealize the transition layer 
to being a thin shell and apply a version of the Sen-Lanczos-Israel junction 
condition formalism [11]. In fact, the latter approaches have been extensively 
analysed in the literature, and applied to a wide variety of scenarios [121 113]. Several 
criteria related to potential observability have been explored [H]. 

The key point of the present paper is to develop an extremely general and 
robust framework that can quickly be adapted to wide classes of generic thin- 
shell gravastars. We shall consider standard general relativity, with gravastars that 
are spherically symmetric, with the transition layer confined to a thin shell. The 
bulk spacetimes (interior and exterior) on either side of the transition layer will 
be spherically symmetric and static but otherwise arbitrary. (So the formalism 
is simultaneously capable of dealing with gravastars embedded in Schwarzschild, 
Reissner-Nordstrom, Kottler, or de Sitter spacetimes, or even "stringy" black hole 
spacetimes. Similarly the gravastar interior will be kept as general as possible for as 
long as possible.) The thin shell (transition layer) will be permitted to move freely 
in the bulk spacetimes, permitting a fully dynamic analysis. This will then allow 
us to perform a general stability analysis, where gravastar stability is related to the 
properties of the matter residing in the thin-shell transition layer. 
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Many of the purely technical aspects of the analysis are very similar in 
spirit to that encountered in a companion paper analyzing thin-shell traversable 
wormholes |15] — mathematically there are a few strategic sign flips — but physically 
the current framework is significantly different. Consequently we shall very rapidly 
find our analysis diverging from the traversable wormhole case. Further afield, we 
expect that the mathematical formalism developed herein will also prove useful when 
considering spacetime "voids" (manifolds with boundary) |16j . 

This paper is organized in the following manner: In Section [2] we outline in 
great detail the general formalism of generic dynamic spherically symmetric thin 
shells, and provide a novel approach to the linearized stability analysis around a 
static solution. In Section |3l we provide specific examples and consider a stability 
analysis by applying the generic linearized stability formalism outlined in section [2j 
Finally, in Section HI we shall draw some general conclusions. Throughout this work, 
we adopt the sign conventions of Misner-Thorne- Wheeler [17] , with c = G = 1. 



2. General formalism 



To set the stage, consider two distinct spacetime manifolds, an exterior A^+, and 
an interior A^„, that are eventually to be joined together across some surface layer 



S. Let the two bulk spacetimes have metrics given by g 



X': 



and g 



X_ 



in terms 



of independently defined coordinate systems x+ and x^. In particular, consider two 
generic static spherically symmetric spacetimes given by the following line elements: 





\ b±{r±y 


dtl + 






r± 







drl 



rldQl. 



(1) 



We take + to refer to the exterior geometry and — to refer to the interior geometry. 
For simplicity we assume that the exterior geometry is asymptotically fiat and define 

i?_ = max{r : 6+(r) = r}. (2) 

(More generally, if the exterior geometry is not asymptotically fiat, the maximum 
should be taken over whatever black hole horizons are present, but excluding the 
cosmological horizons.) Similarly we assume that the interior geometry is regular at 
the origin (r = 0) and define 

i?+ = min{r : 6_(r) = r}. (3) 

Note the (at first glance) counter- intuitive placement of the ± on the quantities R^-: 
The conventions are chosen so that is the furthest outwards one can extend the 
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interior geometry before hitting a horizon, whereas -R_ is the furthest inwards one 
can extend the exterior geometry before hitting a horizon. 

Since the whole point of a gravastar model is to avoid horizon formation we 
certainly desire < -R+. In particular, if the thin-shell transition layer E is located 
at r = a(r), then to avoid horizon formation we demand 

< a(r) < (4) 

The key issue of central interest in this article is the dynamics of this surface layer. 

2.1. Bulk Einstein equations 

Using the Einstein field equation, Gf^^ = Stt T^^, (with c = G = 1), the (orthonormal) 
stress-energy tensor components in the bulk are given by 

Pr{r)= - [2$'(6 - r) + 6'] , (6) 

p^(r) = - Y^^li-b + 3rb' - 2r)<^' + 2r(6 - r)(<l>')^ + 2r{b - r)<^" + b"r] , (7) 

where the prime denotes a derivative with respect to the radial coordinate. Here p{r) 
is the energy density, Pr{r) is the radial pressure, and pt{r) is the lateral pressure 
measured in the orthogonal direction to the radial direction. The ± subscripts were 
(temporarily) dropped so as not to overload the notation. Note that in obtaining 
the individual field equations (I5l)-(I71), instead of using an orthonormal basis in the 
Einstein field equation one could simply, (because of the diagonal form of the metric 
in the current situation), consider the mixed tensor components, i.e., G^'^ = SttT^'^. 

2.2. Null energy condition 

Consider the null energy condition (NEC): T^yk^k'^ > 0, where T^j, is the stress- 
energy tensor and k^^ any null vector. Then along the radial direction, with 
kf^ = (1, ±1, 0, 0) in the orthonormal frame where T^,> = dia.g[p{r),Pr{r),pt{r),pt{r)], 
we have the particularly simple condition 

Tf,o k^ k" = p{r) + pr{r) = ^"""^^ > 0. (8) 

By hypothesis r > h{r) in both the interior and exterior regions of the gravastar, so 
the radial NEC reduces to $'(r) > 0. The NEC in the transverse direction, p+pt > 0, 
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does not have any direct simple interpretation in terms of the metric components. 
In most gravastar models the NEC is taken to be satisfied, though the status of the 
NEC as fundamental physics is quite dubious [18]. 



2.3. Transition layer 

The interior and exterior manifolds are bounded by isometric hypersurfaces S+ and 
S_, with induced metrics g^j and gjj. By assumption g^j{^) = gljiC) = QijiO^ with 
natural hypersurface coordinates ^* = {T,6,(f)). A single manifold Ai is obtained by 
gluing together Ai^ and at their boundaries. SoA^ =A^+UA^_, with the 
natural identification of the boundaries S = S_|_ = S_. The intrinsic metric on S is 

dsl = -dr^ + a(r)2 {dO^ + sin^ 6 d(f)^). (9) 

The position of the junction surface is given by x^{T,6,(j)) = (t{T),a{T),9,(j)), and 
the respective 4-velocities (as measured in the static coordinate systems on the two 
sides of the junction) are 



U^=\ , a,0,Ol. (10) 

a 

The overdot denotes a derivative with respect to r, the proper time of an observer 
comoving with the junction surface. The Israel formalism requires that the normals 
point from to [TT]. The unit normals to the junction surface are 



a 

In view of the spherical symmetry these results can easily be deduced from the 
contractions U^U^ = —1, f/'^n^ = 0, and n'^n^ = +1. The extrinsic curvature, or the 
second fundamental form, is defined as Kij = nfj^.^e^^-^e^f^-y Differentiating n^j^e'^^-^ = 
with respect to , we have 

so that general the extrinsic curvature is given by 
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For a thin shell Kij is not continuous across S. For notational convenience, the 
discontinuity in the second fundamental form is defined as Kij = Kf- — K~-. The 
non-trivial components of the extrinsic curvature can easily be computed to be 



K\^ = -dl-^ + a\ (14) 
a V a 



b± (a)—b'j^ (a)a 



= { '^""^ + ^'^{a)Jl - ^ + ^ , (15) 
{.JV^^^f^^ V a J 

where the prime now denotes a derivative with respect to the coordinate a. 

• Note that is independent of the quantities (^±. This is most easily verified 
by noting that in terms of the normal distance I to the shell S the extrinsic 
curvature can be written as K^j = \d^gij = \n'^d^gij = ^n^drgij, where the last 
step relies on the fact that the bulk spacetimes are static. Then since ggo = r^, 
differentiating and setting r — )■ a we have Kqq = a n^. Thus 

i^V = ^, (16) 

(Jj 

which is a particularly simple formula in terms of the radial component of the 
normal vector, and which easily lets us verify ( lT4l) . 

• For Krr there is an argument (easily extendable to the present context) in 
reference [19] (see especially pages 181-183) to the effect that 

K'^^ = g± = (magnitude of the physical 4-acceleration of the transition layer) . 

(17) 

This gives a clear physical interpretation to K'^^ and rapidly allows one to verify 

m- 

• There is also an important differential relationship between these extrinsic 
curvature components 

^{ae*±irV}=e*±^"r^«- (18) 
dr 

The most direct way to verify this is to simply differentiate, using ( !T^ and 
(|T5l) above. Geometrically, the existence of these relations between the extrinsic 
curvature components is ultimately due to the fact that the bulk spacetimes 
have been chosen to be static. By noting that 
1 n . , ^ . 
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we can also write this differential relation as 



da 



(20) 



2.4- Lanczos equations: Surface stress- energy 

The Lanczos equations follow from the Einstein equations applied to the hypersurface 
joining the bulk spacetimes, and are given by 



{21) 



Here 5^ is the surface stress-energy tensor on S. In particular, because of spherical 
symmetry considerable simplifications occur, namely = diag (k'^^, K^g, k^^,) . The 
surface stress-energy tensor may be written in terms of the surface energy density, 
a, and the surface pressure, V, as S*,- = diag(— cr, V, V). The Lanczos equations then 
reduce to 



An Stt 
From equations f lT^ -f lT5]) . we see that: 

1 



a + 2V = -^. 

47r 



a 



V 



Ana 



Svra 



1 



1 + a? + ad — 



b+{a) 



d^-\ll 



6_(a) 



+ a 



o2 



(22) 



(23) 



b-f-{a)+ab',{a) 



2a 



b+{a) 



+ 



1 + + aa 



+ d2 

b-{a)+ah'_{a) 
2a 



6+ (a) 



+ a? a$' (a) 



-I 6_ (a) I • o 
1 — + 



6_(a) 



+ a? a$'_(a) 



(24) 



and finally 



a + 2P = M = ^ 

An An 



a + 



b-\-{a)~ab'^{a) 



2a2 



d + 



- ^ + d2 

a 

b.(a)-ab'^{a) 
2^ 



+ A 1 



6+ (a) 



d2 $' fa) 



-I b_(a) . .9 

1 — + a^ 



6, (a) 



(25) 



Note that a + 2P has a particularly simple physical interpretation in terms of [g\ , 
the discontinuity in 4-acceleration. (This is ultimately related to the fact that the 
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quantity cr + 2V for a thin shell has properties remarkably similar to those of the 
quantity p + 3p for a bulk spacetime.) Furthermore the surface energy density a 
is independent of the quantities The surface mass of the thin shell is given by 

Independent of the state of motion of the thin shell, we have cr(a) > whenever 
6+ (a) > 6- (a), and cr{a) < whenever 6+ (a) < 6_(a). The situation where a = 
corresponds to 6_(a) = 6+ (a), and is precisely the case where all the discontinuities 
are concentrated in i^T^ while K^g is continuous. This phenomenon, the vanishing of 
cr at certain specific shell radii given by 6_(a) = &+(a), is generic to gravastars but 
(because of a few key sign flips) cannot occur for the thin-shell traversable wormholes 
considered in [151 IlHl ED]- This is perhaps the most obvious of many properties 
differentiating gravastars from the thin-shell traversable wormhole case, although 
wormhole geometries surrounded by thin shells, similar to the cases explored in this 
work, have also been analyzed in the literature [2T1 [22]. 



2.5. Static gravastars 

Assume, for the sake of discussion, a static solution at some ao € (i?_,i?+). Then 

1 



a{ao) = - 



47rao 



6+(ao) 



- Wl- 



57rao 



ao 

b+{ao)+aQb'^(a(,) 
2ao 



Moo) 
ao 



(26) 



1 



+ ,/l_M^ao$V(ao) 



ao 



ao 



1 - 



b_(ao)+aofe'_(ao) 
2ao 

_ b-(ao) 
ao 



b-{ao) 



ao 



ao*-(ao) 



(27) 



and finally 



a(ao) + 2P(ao) 



6+ (a)— afe^ 




2a2 




-(a) 




a 


b-{a)-ah[ 




2a? 




-(a) 
a 



- A 1 ^ a 



(28) 
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(See also equation (32) of [23].) Now taking — )■ R-, we have 

1 



+ 



However taking ag — ?■ -R+ we have 

1 



That is: 



a(/?4 




> 0. 



< 0. 



(29) 



(30) 



(31) 



Applying the mean value theorem, for gravastars there will always be some Rq G 
(-R_,-R+), possibly many such Rq, such that cr{Ro) = 0. (This phenomenon cannot 
occur for the thin-shell traversable wormholes considered in p3| IT9l [20] . because of 
key sign flips — for thin-shell traversable wormholes we always have a < 0.) The Rq 
such that cr(i?o) = is clearly a special place for gravastars. Explicitly this occurs 
when 



b+{Ro) = b-{Ro), 



and in fact for 



< 6+(i?o) = b^Ro) < R+. 



(32) 



(33) 



At this special point the discontinuities are concentrated in K'^^ while K^q is 
continuous. That is 



r{Ro) 



so that 



1 



WnR 







b'_{Ro)-b'4Ro 



b±{Ro) 



1 - 



+ 2W 1 - RolKiRo) - $'_(^o)] 



(34) 



V{Ro 



[K^r] [go] 



Sir Sir 

Even though a = 0, one needs V ^ because of the non-zero 4- acceleration Qq. 

If one also demands that the surface pressure at Rq also be zero, V{Rq) 
one must impose the additional condition 

6V -y_ = 2 {Ro - 6±) - , 



(35) 



(36) 
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which is equivalent to 

(e2*+[l - b^/ao])' = {e'^'ll - bjao])' . (37) 

Ro Ro 

If $'_t(-Ro), then not only does 6+(-Ro) = b-{Ro) but also b'^{Ro) = 6'_(_Ro) — these 
two conditions give a zero pressure and zero density shell. More specifically, in terms 
of standard nomenclature, if the surface stress-energy terms are zero, the junction 
is denoted as a boundary surface; if surface stress terms are present, the junction is 
called a thin shell. 

2.6. Surface stress estimates 

It is interesting to obtain some estimates of the surface stresses. For this purpose, 
consider for simplicity an exterior geometry that is Schwarzschild-de Sitter spacetime, 
so that 

6+(r) = 2M+^r^ <l>+(r) = 0. (38) 

An important quantity that will be play a fundamental role throughout this paper 
is the surface mass of the thin shell, which is given by = Ana'^a. For the exterior 
Schwarzschild-de Sitter spacetime, and considering an arbitrary interior geometry, 
the surface mass of the static thin shell is given by 




m.(ao) = ao h/1 - " V ^ " ^ " I • (39) 

Note that one may interpret M as the total mass of the system, as measured in an 
asymptotic region of the spacetime. Solving for M, we have 





M = ^ + ™,(a.) I ,/ 1 - - ^ j - -a^ (40) 

For the Schwarzschild-de Sitter spacetime A > 0. For the range < 9AM^ < 1, the 

factor = —gu = (1 — 2M/r — Ar^/3) possesses two positive real roots, and r^, 
corresponding to the black hole and the cosmological event horizons: 

n = 2A-1/2 cos(a/3) , (41) 

rc = 2A"^/2 cos(a/3 + 47r/3). (42) 

Here cos a = — 3MA^/^, with vr < a < 37r/2. In this domain we have 2M < < 3M 
and > 3M 
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Defining suitable dimensionless parameters, equations ( I26|) - fl271) take, for tlie 
Scfiwarzscliild-de Sitter solution, the form 




/i(ao) = x { \/l-x b{x) - \/l-x- j , (43) 

1-f- l-'^[xbix) + b'_] 



n(ao) = X , ' ) ' ' - C- - ^^(^) • (44) 

Here we define: x = 2M/a, (3 = 9AM^, b{x) = 6(a)/(2M), /x = SvrMa, and set 
n = IGnMV. In the analysis that follows we shall assume that M is positive, 
M > 0. 

As a specific application, consider the standard gravastar picture, which consists 
of a Schwarzschild exterior geometry and an interior de Sitter spacetime. Thus the 
surface stresses, equations (H3!) - (H^ . are obtained by setting A = 0, (i.e., (3 = 0), 
while 6_(ao) = a^/R^ and $_(ao) = 0. Here we have defined E? = 3/A_, for 
simplicity. For this case the total mass of the system M is given by 




M=^ + m.(ao)|^yi-^-;^l . (45) 

Using the dimensionless parameters, x = 2M/aQ and y = 2M/R, while /i(ao) 
87iMa{ao) and H(ao) = 16nMV{ao), equations (|13D-(144]), take the form 




/i(ao) = -X I y/l-x - Y 1 - — I , (46) 
1-^ l-2< 



H(ao)=x| ^=^-^^^ I . (47) 



The qualitative behaviour is depicted in figure [TJ The left plot represents the surface 
energy density, a. Note that a is positive in the range 2M < < (2Mi?^)^/^, and 
negative for (2Mi?^)^/^ < < R. The right plot represents the surface pressure, 
V. The latter diverges as approaches 2M {R) from the right (left). A detailed 
stability analysis of the thin shell for this specific configuration will be presented in 
section 13.11 
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Figure 1. Schwarzschild-de Sitter gravastar: These plots depict the quahtative 
behaviour of the surface stresses of a gravastar with interior de Sitter and exterior 
Schwarzschild spacetimes. We have considered the dimensionless parameters, 
X = 2M/ao and y = 2M/R, and /i(ao) = SnMaiao) and n(ao) = 16TrMV{ao). 
The left plot represents the surface energy density, a. Note that a is positive in 
the range 2M < oo < {2MR'^y^^, and negative in the range (2Afi?2)V3 < ao < R. 
The right plot represents the surface pressure, V. The latter surface pressure 
diverges as oq approaches 2M (i?) from the right (left). See the text for more 
details. 



2.7. Conservation identity 

The first contracted Gauss-Codazzi equatiorjil is 

G,, n^n' = ^ {K' - K,,K'^ - 'R) . (48) 
The second contracted Gauss-Codazzi equatior(§| is 

G.^ef^^n'^ = - K,, . (49) 
Together with the Lanczos equations this provides the conservation identity 

^ j\i - 

I Sometimes called simply the Gauss equation, or in general relativity more often referred to as 
the "Hamiltonian constraint" . 

§ Sometimes called simply the Codazzi or the Codazzi-Mainardi equation, or in general relativity 
more often referred to as the "ADM constraint" or "momentum constraint" . 



(50) 
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where the convention [X]^ = X+l^ — X is used. When interpreting this 
conservation identity, consider first the momentum flux deflned by 



T 



n 



a\ I - 



b{a) 



+ d2 



1 _ M 



(51) 



where Tjj and Tfr are the bulk stress-energy tensor components given in an 
orthonormal basis. This flux term corresponds to the net discontinuity in the (bulk) 
momentum flux = T^^, which impinges on the shell. Applying the (bulk) 
Einstein equations we see 



T 



n 



Ana 



$'+(a)Wl 



M^ + a^-$'_(a)Jl-^ + d2 



It is useful to deflne the quantity 
1 



+ a 



hi 



(52) 



(53) 



and to let A = Ana? be the surface area of the thin shell. Then in the general case, 
the conservation identity provides the following relationship 



da , 1 dA _ . 



or equivalently 



d{oA) 



V 



dA 



SAd. 



(54) 



(55) 



dr dr 

The flrst term represents the variation of the internal energy of the shell, the second 
term is the work done by the shell's internal force, and the third term represents the 
work done by the external forces. Once could also brute force verify this equation 
by explicitly differentiating ( 123|) using and the relations (fT8|) . If we assume that 
the equations of motion can be integrated to determine the surface energy density 
as a function of radius a, that is, assuming the existence of a suitable function a{a), 
then the conservation equation can be written as 

a' = --{a + V) + E, (56) 

a 

where a' = da /da. We shall carefully analyze the integrability conditions for a{a) 
in the next sub-section. For now, note that the flux term (external force term) S is 
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automatically zero whenever ^± = 0; this is actually a quite common occurrence, for 
instance in either Schwarzschild or Reissner-Nordstrom geometries, or more generally 
whenever p + = 0, so it is very easy for one to be mislead by those special 
cases. In particular, in situations of vanishing flux S = one obtains the so-called 
"transparency condition", [G^^ n"]'^ = 0, see [2l]. The conservation identity, 
equation (l50l) . then reduces to the simple relationship a = —2 [a + V)a/a. But in 
general the "transparency condition" does not hold, and one needs the full version 
of the conservation equation as given in equation fl55|l . 

2.8. Integrability of the surface energy density 

When does it make sense to assert the existence of a function a (a)? Let us start 
with the situation in the absence of external forces (we will rapidly generalize this) 
where the conservation equation, 

& = -2{a + V)d/a, (57) 

can easily be rearranged to 

2-. (58) 



a + V a 

Assuming a barotropic equation of state V{a) for the matter in the gravastar 
transition layer, this can be integrated to yield 

r da r da , , , , 

This implies that a can be given as some function a{a) of a, and by the inverse 
function theorem implies over suitable domains the existence of a function a{a). 
Now this barotropic equation of state is a rather strong assumption, albeit one that 
is very often implicitly made when dealing with thin-shell gravastar (or thin-shell 
wormholes, or other thin-shell objects). As a first generalization, consider what 
happens if the surface pressure generalized is to be of the form V{a, cr), which is not 
barotropic. Then the conservation equation can be rearranged to be 

= J±^Ih^, (60) 

a 

This is a first-order (albeit nonlinear and non-autonomous) ordinary differential 
equation, which at least locally will have solutions cr(a). There is no particular 



Generic thin-shell gravastars 



16 



reason to be concerned about the question of global solutions to this ODE, since in 
applications one is most typically dealing with linearization around a static solution. 

If we now switch on external forces, one way of guaranteeing integr ability 
would be to demand that the external forces are of the form E{a,a), since then 
the conservation equation would read 



which is again a first-order albeit nonlinear and non-autonomous ordinary differential 
equation. But how general is this S = S(a, a) assumption? There are at least two 
nontrivial situations where this definitely holds: 

• If $+(a) = $-(a) = then S = — $'(a) a, which is explicitly of the required 
form. 

• If 6+ (a) = 6_(a) = b{a), but the ^± are unequal, then a = regardless of the 
location and state of motion of the transition layer. Furthermore S = 2V/a. 
(This would make for a somewhat unusual gravastar.) 

• If both 6+(a) = 6_(a) = b{a) and $+(a) = ^-(a) = $(a), then the situation is 
vacuous. There is then no discontinuity in extrinsic curvatures and the thin shell 
carries no stress-energy; so this in fact corresponds to a "continuum" gravastar 
with b{r)/r < 1 for all r e (0, oo). 

But in general we will need a more complicated set of assumptions to assure 
integrability, and the consequent existence of a function a{a). A model that is always 
sufficient (not necessary) to guarantee integrability is to view the exotic material in 
the transition layer as a two- fluid system, characterized by a± and V±, with two 
(possibly independent) equations of state V±{a±). Specifically, take 




(61) 



a 




(62) 



(63) 



In view of the differential identities 



(64) 



each of these two fluids is independently subject to 



(65) 
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which is equivalent to 



2e*± 



(66) 



With two equations of state V±{a±) these are two nonhnear first-order ordinary 
differential equations for a±. These equations are integrable, implicitly defining 
functions (T±{a), at least locally. Once this is done we define 

0"(a) = cr+(a) — cr_(a), (67) 

and 

ms{a) = 4:7ia{a) a^. (68) 

While the argument is more complicated than one might have expected, the end 
result is easy to interpret: We can simply choose cr(a), or equivalently ms{a), as an 
arbitrarily specifiable function that encodes the (otherwise unknown) physics of the 
specific form of matter residing on the gravastar transition layer. 



2. 9. Equation of motion 

To qualitatively analyze the stability of the gravastar, assuming integrability of the 
surface energy density, (that is, the existence of a function o"(a)), it is useful to 
rearrange equation (!23l) into the form 

-a^ + V{a)=0, (69) 



where the potential V{a) is given bjl^l 



V{a) 



1 



b{a) 



2a 



1 2 



AW 

ms{a) 



1 2 



(70) 



Here ms{a) = Arra^ a{a) is the mass of the thin shell. The quantities b{a) and A(a) 

are defined, for simplicity, as 

- &+(a) + &_(a) 

= n ' i^^) 



A(a) 



2 

fa) - 6_ 



(72) 



f This equation only valid for cr ^ due to a divide-by-zero problem. The a = case is a special 
one worth separate consideration: 

2~1 



1 



1 



bia) 



8ttV 
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respectively. This gives the potential V{a) as a function of the surface mass ms{a). 
By differentiating with respect to a, (using (ITQll ). we see that the equation of motion 
implies 

a = -V'{a). (73) 

It is sometimes useful to reverse the logic flow and determine the surface mass as a 
function of the potential. Following the techniques used in [101 US], suitably modified 
for the present context, a brief calculation yields 



6+ (a) 



2V{a) 



bja) 



2V{a) 



(74) 



with the negative root now being necessary for compatibility with the Lanczos 
equations. Note the logic here — assuming integrability of the surface energy density, 
if we want a specific V{a) this tells us how much surface mass we need to put on 
the transition layer (as a function of a), which is implicitly making demands on 
the equation of state of the matter residing on the transition layer. In a completely 
analogous manner, the assumption of integrability of a{a) implies that after imposing 
the equation of motion for the shell one has 



(j(a) 



Area 



6+ (a) 



2V{a) 



(75) 



while 



V 



Sna 



1 - 2V{a) - aV'{a) 



fe+ {a)-\~ab'^ (a) 
2a 



1 - 



6+ (a) 



1 - 2V{a) - aV'{a 



2V{a) 

b^{a)+ab'_{a) 



+ \ I 



6+ (a) 



2V{a) a$V(a) 



2a 



b-{a) 
a 



2V{a) 



b^(a) 



2V{a) a^'_{a) 



and 



4:7ia 



l-M^-2na)-$'_(a)A/l-^ 



2V{a) 



(76) 



(77) 



The three quantities {a{a),V{a),E{a)} (or equivalently {ms{a),V{a),E{a)}) are 
related by the differential conservation law, so at most two of them are functionally 
independent. 
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2.10. Linearized equation of motion 

Consider a linearization around an assumed static solution (at ao) to the equation 
of motion + V{a) = 0, and so also a solution of a = —V'{a). Generally a Taylor 
expansion of V{a) around gq to second order yields 

V{a) = V{ao) + V'{ao){a - oq) + ^V"{ao){a - aof + 0[{a - aof] . (78) 

But since we are expanding around a static solution qq — do — 0, we automatically 
have V{ao) — V'{ao) — 0, so it is sufficient to consider 

V{a) = ^V"{ao){a - aof + 0[{a - aof] . (79) 

The assumed static solution at is stable if and only if V{a) has a local minimum 
at Oq, which requires V"{aQ) > 0. This will be our primary criterion for gravastar 
stability, though it will be useful to rephrase it in terms of more basic quantities. 

For instance, it is extremely useful to express mg{a) and m"{a) by the following 
expressions: 

ms{a) a j JbM/fl+^n±_ (fe-(«)/«)^ + 2^ (a) , 

mJa) — H \- t: S — = = r : (oUj 

^ a 2 \^l-b+{a)/a-2V{a) ^Jl - b_{a)/a - 2V{a) ' ^ ^ 

and 

, ; (bJa)/ay + 2V'{a) ib4a)/ay + 2V'{a] 



^yl-b+{a)/a-2V{a) ^Jl ~ b_{a)/a - 2V{a) 
a\ \{b^{a)la)' ^2V\a)f \{b^{a)la)' ^2V\a)f \ 
4 \ [1 - 6+(a)/a - 2V{a)fl'^ [1 - 6_(a)/a - 2\/(a)]3/2 / 

l«[ {hMla)" + 2V'\a) _ {b {a)lar^2V'\a) \ 
2 \yi-6+(a)/a-2y(a) ^1 - 6_(a)/a - 2\/(a) J ' 

Doing so allows us to easily study linearized stability, and to develop a simple 
inequality on m'/(ao) by using the constraint K"(ao) > 0. Similar formulae hold 
for a' (a), o""(a), for 'P'(a), V"{a), and for S'(a), S"(a). In view of the redundancies 
coming from the relations ms{a) = 4na{a)a'^ and the differential conservation law, 
the only interesting quantities are S'(a), S"(a). 
It is similarly useful to consider 



47rS(a) a 



^'^{a)\ll- ^ - 2V{a) - $'_(a)7l - ^ - 2V{a) 



(82) 
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for which an easy computation yields: 

[47rS(a)a]'= + |$'|(a)Vl - b+{a)/a - 2V{a) - ^'L{a)^/l - b^{a)/a - 2F(a)} 



1 I {b+{a)/ay + 2V'{a) _ (6-(a)/a)' + 2V'{a) 



2[^' ' ^l-b+{a)/a-2V{a) \/l - b-{a)/a - 2V{a) 

and 

[47r S(a) a]" = {^'^{a)^/l - b+{a)/a - 2V{a) - ^'!!{a)^/l - b-{a)/a - 2V{a)^ 



(83) 



[(b+(a)/ay + 2y-(«)]^ _ J(M^)/ay + 2n^l 

[1 - b+{a)/a - 2y(a)]3/2 ^ [1 _ 6_(a)/a - 2y(a)]3/2 / 



1 (6+(a)/ar + 2y--(a) _ ^, (6_(a)/a)" + 2y"(a) 



2 [ ' ^/l-b+{a)/a-2V{a) - b-{a)/a - 2V{a) 

We shall now evaluate these quantities at the assumed stable solution qq- 
2.11. The master equations 

In view of the above, to have a stable static solution at Oq we must have: 



(84) 



) , (85) 



while 



<iao)-^-li - /-^-(^°) 1. (86) 

2ao 2 ^1 - 6^(ao)/ao a/I - 6_(ao)/ao J 

The inequality one derives for m"{aQ) is now trickier since the relevant expression 
contains two competing terms of opposite sign. Provided 6+(ao) > 6_(ao), which is 
equivalent to demanding cr(ao) > 0, one derives 

1 r [b+{ao) - ao6V(«o)]^ [^-(oo) - aob'_{ao)f 



m':{a,) > + ^ { 



[l-6+(ao)/ao]3/2 [1 - 6_(ao)/ao]3/2 
1 I 6': fan) 6" fan) 
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However if 6+(ao) < 6_(ao) the direction of the inequahty is reversed. This last 
formula in particular translates the stability condition V"{aQ) > into a rather 
explicit and not too complicated inequality on m^'(ao), one that can in particular 
cases be explicitly checked with a minimum of effort. 

In the absence of external forces this inequality is the only stability constraint 
one requires. However, once one has external forces (H which requires $± 0), 
there is additional information: 

[47rH(a)a]'L^ = + {$^(a) ^1 - b+{a)/a - <^'L{aWl - &_(«)/«} 



^yl - b+{a)/a 



^\ — b_{a)/a 



(88) 



ao 



Provided (oq) / ^/^ — b+{ao)/ao > $'_(ao)/-\/l — 6_(ao)/ao, we have 
[47rS(a)a]"|„„ < ^1 - - ^-(«) Vl - &-(«)/«} 

(6+(a)/a)' ^ (6_(a)/a)' 




1 — b^{a)/a 
[{b4a)/ayr 



[1 — b^{a)/a 



3/2 



a 



1^1 — 6_ (a) /a 
[(6_(a)/a 




-$'_(a)- 



[l-6_(a)/a]3/2 
(6_(a)/a)" 



ao 



(89) 



ao 



b^{a)/a ' ' -^1 — b^{a)/a 

If $+(ao)/-\/l — 6+(ao)/ao < $'_(ao)/-\/l — 6-(ao)/ao then the direction of the 
inequality is reversed. Note that these last two equations are entirely vacuous in 
the absence of external forces, which is why they have not appeared in the literature 
until now. 



3. Specific gravastar models 

In discussing specific gravastar models one now "merely" needs to apply the general 
formalism described above. Up to this stage we have kept the formalism as general 
as possible with a view to future applications, but we shall now focus on some more 
specific situations. 
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3.1. Schwarzschild exterior, de Sitter interior 



The traditional gravastar, first considered in [T] and [2], has a Schwarzschild exterior 
(with &+(r) = 2M and $+(r) = 0) and a de Sitter interior (with 6_(r) = r^ / B? 
and $_(r) = 0), but with a complicated transition layer. Thin-shell Schwarzschild- 
de Sitter gravastars were first explicitly discussed in [10]. The parameters are chosen 
such that the transition layer is located at some 2M < a < R. (So the transition 
layer is situated outside the region where the Schwarzschild event horizon would 
normally form, and inside the region where the de Sitter cosmological horizon would 
form). One normally is rather noncommittal regarding the physics of the transition 
region, however, for the present case, the surface stresses of the thin shell are given 
by 



a 



V 



Ana 



2M 



Sua 



1 + + aci - — 



i?2 



a 



1 + + ad - 2^ 



7? 



(90) 



(91) 



The external forces vanish (S = 0), as $± = 0, and a > for a(r) < (2Mi?2)^/^, as 
can be readily verified from equation ( 190|) . To have a stable static solution at ao we 
must have: 



a(ao) 



47rao 



1 



2M 
ao 



57rao 



M 
ao 



1 



o °0 



i?2 



\ _ 2M 

ao 



1 _ 



(92) 



(93) 



with 2M < ao < -R, a situation which has already been extensively analysed in 
Section 12.61 As pointed out in the general discussion, a takes finite values with 
different signs in the endpoints of the range between 2M and i?, being positive for 
2M < ao < {2MB?YI^ and negative for {2MR^Y/^ < < R. The surface pressure 
V tends to +oo when ao approaches 2M (i?) from the right (left). It can be seen 
that V never vanishes in this interval, and that its derivative with respect to ao tends 
to — oo (+oo) when ao goes to 2M (i?) from the right (left). That is, V (ao) evolves 
from infinitely large values, to a minimum non-vanishing value, before then again 
going to infinity. 
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Now consider the mass 

m.(ao) = -ao <! a/1 - — - \/l - ^ ^ (94) 
while 

\ ^1 - 2M/a„ 

The inequahty one derives for m"(ao) is now trickier since the relevant expression 
contains two competing terms of opposite sign. However, considering that the 
physical solution should have cr > (equivalent to ao < (2Mi?^)^/^), one finds 
that stability requires 

We can recast this as 

(M/o^y' (a„/-R)' [3-2(a„/fl)'] 
<•» ^ 11 - 2M/ao]3/^ 11 - ia^lKffl- ' 

In figure [2] we show the surface which is produced when this inequality saturates, 
where the stability region are represented above this surface. It is interesting to note 
that the possible positions of a static thin shell were studied in reference [10]. As 
those solutions, ao, were obtained by considering a particular equation of state for 
the matter on the shell, and a and V are independent of V" ^ then the solutions would 
be the same for the zero potential as for the linearized potential. The only difference 
between the two models is that whereas if \^ = then the solution would be on the 
surface depicted in figure Ej if we consider the linearized potential then the solutions 
should be in the region above the surface in order to assure stability of the static 
solution (which is equivalent to demanding V" (ao) > 0). 



3.2. Bounded excursion gravastars 

On the other hand, once one obtains a static solution at ao by requiring a specific 
equation of state (for example, stiff matter on the shell a = V), and verifies 
the stability by checking that the second derivative of the mass of the thin shell 
is in the stability region, it is easy to obtain dynamic "stable" solutions of the 
'bounded excursion' type by deforming the linearized potential. Thus, considering 
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Figure 2. Schwarzschild-de Sitter gravastar: The function ao m" in the case that 
= 0. The surface is only defined for values of Oq < i? and 2M < oq, as expected. 
Models producing a function ao m" above this surface would be in the stability 
region. 



V{a) = ^ (a — ao) — with e sufficiently small, one can obtain the equation of 
motion of the shell; this is 

a(r) = ao + - sin [7 (r - To)] . (98) 

7 

Therefore, the shell expands from a minimum size with ai = ao — e/7 to a maximum 
size corresponding to a2 = ao + e/7. It then contracts to ai, starting a new cycle of 
evolution after reaching this value. Now, it can be clearly understood what we mean 
with e sufficiently small, because this behavior makes sense for a stable gravastar 
only if ai > 2M and a2 < R, which implies e < 7 (ao — 2M) and e < 7 (i? — ao), 
respectively. In figure [3] we show the behavior of the potential, which is related with 
d through the equation of motion, as a function of a. This potential vanishes at ai 
and a2, where d = but d = ±67, respectively. This acceleration imposes that the 
shell rolls down the potential when it has reached both its minimum and maximum 
sizes. The shell continues evolving when its radius takes the value ao because at this 
point d 7^ 0. 

Note that since 

dr ^ l-b±{a)/a _ 
dt± ^l-6±(a)/a + d2' 
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Figure 3. Bounded excursion potential obtained by deforming the linearized 
potential. V{a) vanishes at ai^2, with ai < oq < 02, being ao the minimum of 
the potential. The shell expands from ai to 02, and then it contracts again due to 
the nonvanishing acceleration at this point. 

one would have, in general, a different equation of motion of the shell in terms 
of the time coordinate of the exterior geometry and the time coordinate of the 
interior geometry, that is a(t+) 7^ a(t_). Nevertheless, as r can be considered as 
a parameter in both regions of the space, the shell radius would always be bounded 
by Oi and 02- Thus, if ai and 02 can be reached at a finite t±, then the shell would 
be vibrating, although with a different kind of vibration as seen using t+ or t_. (See, 
for instance, [25].) 

Finally, we should comment on some features of the material on the shell. 
Although we are deforming a solution corresponding to a particular kind of material 
on the shell, that is with a given equation of state parameter relating a and V, 
the corresponding dynamic solution would not have a constant equation of state 
parameter, at least in the general case, because the surface stresses, a{a) and V{a), 
have a different dependence on the trajectory (see equations ( 190|) ). Therefore, the 
shell would be filled by material changing its behavior during the evolution of the 
shell. Moreover, even if the original static solution corresponds to a material on the 
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shell fulfilling the energy conditions, its dynamic generalization could easily violate 
those conditions at some stage of the evolution of the shell. Thus, one should carefully 
study that this is not the case in each particular model. In figure IHwe have depicted 
the equation of state parameter {w = V/a) a.s a function of time, and the energy 
density on the shell as a function of the pressure for a particular model obtained by 
deforming a stable static solution with dust matter on the shell. It can be seen that 
the pressure on the dynamic shell never vanishes, although the solution is obtained 
by deforming one with 'P(ao) = 0. 



w 




Figure 4. Bounded excursion: We consider the deformation e ~ 0.01 of a stable 
static solution with dust matter on the shell, M = 1, oq = 4.178821374980832, 
R = 22.3607 and 7 = 1. In the left figure we show the evolution of w in terms of 
T. It decreases from its maximum value to its minimum value from ri — T{a = ai) 
to T2 = T{a = 02) and then increases to its maximum value for the next cycle. 
Evolution of the energy density in terms of the pressure is depicted in the right 
figure. Consideration of more than one cycle would lead to the same graphic. 



3.3. Close to critical: Non-extremal 

Another common feature of gravastars is that the transition layer is typically taken 
to be "close" to horizon formation. That is b±{a)/a ~ 1. In fact, it is useful to take 
a linear approximation 

b±{a)/a = 1 T 7±(a - R^) + 0{[a - R±\^) , (100) 

where R- is where the "black hole" horizon would have formed in the exterior 
spacetime, R^ is where the "cosmological" horizon would have formed in the interior 
spacetime. We take the transition layer to be at a position a such that R- < a < R+, 
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and 7± > because we are (for now) avoiding "extremal" geometries. Dismissing 
second order terms and considering vanishing external forces ($±(r) = 0), we have 
1 



cr 



V 



Ana . 



^7+ [a - 



a) + 



+ aa + 7+ (3a/2 - R_) + ad + 7_ (i?+ - 3a/2) 



(101) 
(102) 



(a - R-) + d2 y/7_ (i?+ - a) + d2 

Therefore, cr(a) > for a(r) < (7+-R_ + 7_-R+) /(7+ + 7-). 

This "close to the horizon" approximation would be valid only if the trajectory 
of the transition layer, which can be obtained from the equation of motion (or 
equivalently considering some equation of state), is always in the region i?_ < a < 

In order to analyze the accuracy of this approximation, we consider a stable 
static solution, implying 

1 



cr(ao) 
^(«o) 



47ran 



a/7+ (ao - R-) - \/l-{,R+ - ao) 



7+ (3ao/2 - R_) 7_ (i?+ - 3ao/2) 



V 1+ («o 



Consider stiff matter on the shell, a (oq) 



R.) V7_ {R+ - ao) 
= "P (flo)- Therefore, we have 



7^^ (3/?+ - 7ao/2) 7!^" (7ao/2 - 3i?, 



1/2 



Squaring both sides, and defining the dimensionless quantities a = ao/R. 

and r_ = 7_-R_, with a > 1 and < /3 < 1, one has 



R-/R+, 



49 



+ 



21 

J 



49 

T 



9 



/3 



+ 21 



9 /r 



a — 



/3 V /3 



(103) 
(104) 

(105) 

/3 = 

(106) 



0. 



Thus, by considering the expansion of the background geometries close to where the 
horizon would be formed, we have reduced the problem of finding static and stable 
solutions (which usually involves some highly nontrivial equation) to solving a cubic, 
which can be done analytically. 

Nevertheless, some comments are in order. In the first place, as the RHS of 
equation f llOSp is always positive, the solutions of equation ( I106p would correspond 
to solutions of our problem only if i?+ > 7ao/6, that is /5 < 6/7. In the second 
place, one can consider that the approximation would not be accurate enough if 



Generic thin-shell gravastars 



28 



(ao — -R-)/-R- < 1 and (-R+ — ao)/ao < 1 are not satisfied; thus, the solutions 
would be reliable only if a < 2 and 1/4 < /3. In summary, we should consider 
l/4</3<6/7to solve equation (11061) . and give physical meaning only to solutions 
with 1 < a < 2, if any. In fact, one can see that the results that can be obtained by 
using this approximation are compatible with those of a stiff matter gravastar with 
a Schwarzschild exterior and a de Sitter interior when 1 < a < 2 (see equation (60) 
of reference [10] for the equation that must be solved in that case). 

On the other hand, we can study the stability of the static solutions. In this 
case the inequality ( IHTl) . for a > 0, leads to 

7i (3ao/4 - R^) 7! {R+ - 3ao/4) 
(ao) > TjK ■ (107) 

bAao-R-)f' [7-(i?+-ao)r 
It is more useful to consider m'/ (oq) /7+, which can be written in terms of the 
dimensionless quantities previously introduced, and is given by 



m'^jao) ^ 3a/4-l /r_ 1//3 - 3a/4 
7+ " (a-lf^ Y r+(i//5_a)3/2- 

In figure Owe have drawn this function for the case that the inequality saturates for 
two particular values of ao- Thus the region above the surface corresponds to the 
stability region. 




Figure 5. Close to critical (non-extremal): The left and right figures correspond 
to Og = l.li?_ and Og = 1.5i?_, respectively. It can be noticed that the function 
is not defined for values of R+ < ao, since we must consider R- < ao < R+. 
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3.4- Close to critical: Extremal 



We now consider a gravastar with a transition layer "close" to where the horizon 
would be expected to form, when both bulk geometries are "extremal" . This implies 
that the coefficient of the ffist order term in the expansion f llOOp vanishes. Then, 
assuming that the coefficient of the second order term is not vanishing, one has 

b±{a)/a = 1 - 7|(a - R^f + 0{[a - R±f) , (109) 

with R_ < a < and 7± > 0. Note particularly that we now have a different 
definition for 7-1- . 

Again considering the situation $±(r) = 0, we can express the quantities related 
to the thin-shell as 
1 



a 



V 



Ana 



72 (a-R_y + a^-JjUR- 



a)2 + d2 



(110) 



Sua 



+ ad + 7^ (2a^ + i?^ - 3R^a) d'^ + ad + 7^ {2a'^ + i?+ - 3R+a) 



72 (a - R-f + a2 



a)^ + 



(111) 

Thus, we recover that a(r) < (7+-R_ + 7--R+) /(7+ + 7-) for a{a) > 0, although 7± 
would be different than in the former case analyzed. Considering a static solution, 
one obtains 

a (ao) ^ - 73- [7+ («o - R^) - 1-{R+ - ^o)] , fll2) 



AiraQ 



V{ao) 



7+ {2al + Rl- 3i?_ao) 7_ {2al + Rl - 3R+ao) 



(113) 



Oq — i?_ — Oq 

where we have taken into account R- < a < R^, when simplifying. Following a 
similar procedure to that considered in the non-extremal case, one can see that a 
stiff matter shell must fulfil the equation 

7- (4ag + 3Rl - 7R+ao) ^ 7+ (4ag + 3Rl - 7i?_ao) ^^^^^ 

— do — 

which, taking the same definition of F-t, (3 and a introduced in the former case, leads 
to 



4 (F+ + F_) 



4(F_ + ^)+7(^- + F 



F_) + 3( ^ + F, 



/3 



(3 



a 



a 



(115) 



3 /F 



0. 
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This is again a cubic equation. However, whereas for non-extremal geometries 
we have obtained equation (11061) by squaring equation (11051) . in this case that 
step was not necessary. Thus, all solutions of equation (11151) are also solutions of 
equation (I114p . Therefore, noticing the values of /3 and a for which the approximation 
would be sufficiently accurate, we should solve equation (I115P for 1/4 < /3 < 1 and 
then consider only the solutions with 1 < a < 2. 

A special characteristic of these static solutions can be seen when studying their 
stability. Considering a > 0, the inequality flHTj) can be simplified to obtain 



^>-2fi+^y (116) 

7+ V 7+/ 

which is independent of Oq. As the position ao can be obtained by considering 
a particular equation of state for the material on the shell, this implies that the 
stability of the solution is independent of that equation of stat^. In figure |6l we 
show the behavior of this function. 



2 4 6 8 10 7+ 




Figure 6. Close to critical (extremal): We show the function m"7+ in the case 
that V — 0. The stability region would be above this curve, independent of the 
value of flQ. 

Finally, some comments about the case with one non-extremal and one extremal 
background geometries are in order. In this case, for a stiff matter gravastar, 
one would obtain an equation with the RHS of equation (I105p and the LHS of 

f This interpretation could be reinforced by noticing that in this case m"/7-|_ depends only on the 
geometry and on the potential, which, at the end of the day, implies a dependence on the geometry 
and on a (not onV), ii a ^ 0. 
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equation flll4p . or vice versa. Thus, after squaring, the analogue of equation fll06p 
or f lllSp would be a quintic. Therefore, the consideration of the close-to-horizon 
approximation in this situation would generally not significantly simplify the problem 
of obtaining static solutions, because we are only reducing the polynomial equation 
by one degree (see also reference [10]). 



3.5. Charged dilatonic exterior, de Sitter interior 

An interesting solution to consider is that of an interior de Sitter spacetime, (where 
the metric functions are given by 6_(r) = r^ /B? and $_(r) = 0), while the exterior 
spacetime is given by the dilaton black hole solution, which corresponds to an 
electric monopole. The latter de Sitter-charged dilatonic gravastar is considered here 
for the first time. In Schwarzschild coordinates, this exterior solution is described 
by [261 [271128], 

dsl= - 1 , ] dt^+il , ^ dr^ 

+ r\de^ + sin^ e dip^) . (117) 

Here we have dropped the subscripts -|- for notational convenience. The Lagrangian 
that describes this combined gravitational-electromagnetism-dilaton system is given 
by [261 [271 128] 

C = ^{-R/8n + 2(V^)2 + e-2V'FV47r} . (118) 

(Note that the first charged dilatonic solutions, including black hole solutions, were 
considered by Bronnikov et al. [29]). The non-zero component of the electromagnetic 
tensor is given by F^f, = Q/r"^, and the dilaton field is given by e^^ = 1 — 
QVM(/3 + ^r2 + /32). The parameter (3 is defined by /3 = Q72M. 

In terms of the formalism developed in this paper, the metric functions of the 
exterior spacetime are 




^ 2M 



(119) 



l3 + ^/r^ + P\ 

<l>+(r) = -^ln(l + ^). (120) 
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An event horizon exists at rj, = 2M^\ — (3/M. Note that $+(r) is given by 

which is positive (and so satisfies the NEC) throughout the spacetime. Thus, taking 
into account that $_ = 0, we verify that the stabihty condition imposed by the 
presence of the fiux term is governed by inequahty ( |89|) . 

The thin shell is placed in the region < a < -R, so that > &_(a) (and 

consequently o" > 0) is satisfied. Consequently the stability regions are governed by 
inequality fl57|) . The expressions for inequalities fl57|) and dHH]) are extremely lengthy, 
so that rather than write them down explicitly, we will analyse the qualitative 
behaviour of the stability regions depicted in figures [3-[9l 

Consider as a first example the case for /3/M = 1/2, so that the stability regions, 
governed by the inequalities f l87|) and f l89l) . are depicted in figure [3 The left plot 
describes the stability regions, governed by f lHTl) . which lie above the surface. The 
right plot describes the stability regions depicted below the surface, provided by 
inequahty ( l89l) . As a second case, consider the value /3/M = 3/4, depicted in figure 
[HI We verify that the qualitative results are similar to the specific case of /3/M = 1/2, 
considered above. That is, the stability condition dictated by both the inequalities 
( IHTl) show that the stability regions decrease significantly as one approaches the black 
hole event horizon. The final stability regions are depicted in figure [H and are shown 
in between the shaded regions. Note that the respective stability regions increase for 
increasing values of /3/M. 



4. Summary and Discussion 

In this work, we have developed an extremely general and robust framework leading 
to the linearized stability analysis of dynamical spherically symmetric thin-shell 
gravastars, a framework that can quickly be adapted to wide classes of generic thin- 
shells. We have built on a companion paper, which analyzed thin-shell traversable 
wormholes [15]. We emphasize that mathematically there are a few strategic sign 
fiips, which physically implies significant changes in the analysis. An important 
difference is the possibility of a vanishing surface energy density at certain specific 
shell radii, for the generic thin-shell gravastars considered in this work. Due to the 
key sign flips, ultimately arising from the definition of the normals on the junction 
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Figure 7. Dilaton-de Sitter gravastar: /3/M — 0.5. The left plot describes the 
stability regions above the surface. The right plot describes the stability regions 
below the surface. See the text for details. 




Figure 8. Dilaton-de Sitter gravastar: /3/M = 0.75. The left plot describes the 
stability regions above the surface. The right plot describes the stability regions 
below the surface. See the text for details. 



interface, the surface energy density is always negative for the thin-shell traversable 
wormholes considered in [15]. We have also explored static gravastar configurations 
to some extent, and considered the generic qualitative behavior of the surface stresses. 
Relative to the conservation law of the surface stresses, we have analysed in great 
detail the most general case, widely ignored in the literature, namely, the presence 
of a flux term, corresponding to the net discontinuity in the (bulk) momentum flux 
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2M 




Figure 9. Dilaton-de Sitter gravastar: The final stability regions for a thin shell 
dilaton gravastar for (3/M — 0.5 and (3/M = 0.75 depicted in the left and right 
plots, respectively. The stability regions are represented above the green surface, 
which is given by inequality (|87|: and below the blue surface, given by inequality 
((89|) . Thus the final stability regions are given in between both surfaces. Note that 
the stability regions increase for increasing values of (3/M. 

which impinges on the shell. Physically, the latter flux term can be interpreted as 
the work done by external forces on the thin shell. 

In the context of the linearized stability analysis we have reversed the logic 
flow typically considered in the literature and introduced a novel approach. More 
specifically, we have considered the surface mass as a function of the potential, so that 
specifying the latter tells us how much surface mass we need to put on the transition 
layer. This procedure implicitly makes demands on the equation of state of the 
matter residing on the transition layer and demonstrates in full generality that the 
stability of the gravastar is equivalent to choosing suitable properties for the material 
residing on the thin shell. We have applied the latter stability formalism to a number 
of specific cases, namely, to the traditional gravastar picture, where the transition 
layer separates an interior de Sitter space and an exterior Schwarschild geometry; 
bounded excursion; close-to-horizon models, in both the extremal and non-extremal 
regimes; and finally, to specific example of an interior de Sitter spacetime matched to 
a charged dilatonic exterior. This latter case is particularly interesting, as it involves 
the flux term. 

In conclusion, by considering the matching of two generic static spherically 
symmetric spacetimes using the cut-and-paste procedure, we have analyzed the 
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stability of thin-shell gravastars. The analysis provides a general and unified 
framework for simultaneously addressing a large number of gravastar models 
scattered throughout the literature. As such we hope it will serve to bring some 
cohesion and focus to what is otherwise a rather disorganized and disparate collection 
of results. 
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